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Abstract 

In some quasi-one-dimensional weakly disordered media, impurities are large and rare 
rather than small and dense. For an Anderson model with a low density of strong im- 
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f-^. ■ purities, a perturbation theory in the impurity density is developed for the Lyapunov 

CN . exponent and the density of states. The Lyapunov exponent grows linearly with the den- 

sity. Anomalies of the Kappus-Wegner type appear for all rational quasi-momenta even 
O ' m lowest order perturbation theory. 
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43 1 Introduction 
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A perturbative formula for the Lyapunov exponent of a one-dimensional random medium for 
weakly coupled disorder was first given by Thouless |Thoj and then proven rigorously by Pastur 
and Figotin jPFj. Anomalies in the perturbation theory at the band center were discovered by 
Kappus and Wegner |KWj and further discussed by various other authors |DG[ IBK| ICKj . The 
Lyapunov exponent is then identified with the inverse localization length of the system. This 
short note concerns the behavior of the Lyapunov exponent for a low density of impurities, 
each of which may, however, be large. The presented method is as |.ISS| ISch[ ISSSj a further 
application of diagonalizing the transfer matrices without perturbation (here the low density 
of impurities) and then rigorously controlling the error terms by means of oscillatory sums 
of rotating modified Priifer phases. Some of the oscillatory sums remain large if the rotation 
phases (here the quasi-momenta) are rational. This leads to supplementary contributions of 
the Kappus-Wegner type. 

The calucalations are carried through explicitely for the one-dimensional Anderson model, 
but the method transposes also to more complicated models with a periodic background as well 
as low-density disorder with correlations similar to the random polymer model |JSS . Extension 
to a quasi-one- dimensional situation as in |Schj should be possible, but is even more cumbersome 
on a calculatory level. It is also straightforward to calculate and control higher order terms in 
the density. 

As one motivation for this study (apart from a mathematical one), let us indicate that 
a low density of strong impurities seems to describe materials like carbon nanotubes more 
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adequately than a small coupling limit of the Anderson model. Indeed, these materials have 
perfect cristaline structure over distances of microns which leads to a ballisistic transport over 
such a distance |FPWD| . The existing few defects are, on the other hand, quite large. Coherent 
transport within a one-particle framework should then be studied by a model similar to the 
one considered here. However, it is possible that the impurties rather play the role of quantum 
dots so that Coulomb blockade is the determining effect for the transport properties [IMBCYL] 
rather than the coherent transport studied here. 



2 Model and preliminaries 

The standard one-dimensional Anderson Hamiltonian is given by 
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Here uo = (u n )nez is a sequence of independent and identically distributed real random variables. 
The model is determined by their probability distribution p depending on a given density 
pe[0,l]: 

p = (l - p)s Q + p p ; (1) 

where p is a fixed compactly supported probability mesure on R. This measure may simply be 
a Dirac peak if there is only one type of impurity, but different from £ . Set £ = supp(p) and 
£ = supp(p). The expectation w.r.t. the p's will be denoted by E, that w.r.t. the p's by E, 
while E,, and E„ is the expectation w.r.t. p and p respectively over one random variable ueE 
only. 

In order to define the Lyapunov exponent, one rewrites the Schrodinger equation H^ip = Eifj 
using transfer matrices 

VVh-i \ _ t e ( \ t e _ ( v n - E -1 



We also write for if v n = v. Then the Lyapunov exponent at energy E 6 M associated 
to products of random matrices chosen independently according to p from the family {Tf) vG Y, 
of SL(2, R)-matrices is given by 

i( P ,E) = ^^^(lln^ll) ■ ( 2 ) 

The aim is to study the asymptotics of 7(p, E) in small p for \E\ < 2. 

In order to state our results, let us introduce, for E = — 2 cos(A;) with k e (0,7r), the basis 
change M e SL(2,M) and the rotation matrix Rk by the quasi-momentum k: 

If sm(k) 0\ ( cos(A;) — sin(fc) 

M = — . . /,\ , , Rk " 



y/sm(k) V _ cos ( k ) 1 / ' ' " V sin ( k ) cos ( /c ) 
It is then a matter of computation to verify 
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Next we introduce another auxiliary family of random matrices. Set E = [— f , f ) x E and, 
for (if), v) G E: 

f E v = /.',, .1/7/" .1/ ; . 

The following probability measures on E will play a role in the sequel: p^ = ® p and 
p tf = ( - Y^-i 5zl ) <8> p for o G N. The Lyapunov exponents associated to these families 

\1 '—'V— 1 l\q Iql J 

of random matrices are denoted by 7oo (E) and 'jq(E) respectively. It is elementary to check 
that the subgroups generated by matrices corresponding to the supports of poo and p q are non- 
compact, which implies |BLj that the corresponding Lyapunov exponents are strictly positive. 

The matrices T E and T E ^ V induce actions Se, v and Se,^, v ° n ^ via 

MTfM-'ee f E v e e 
es BtV w \\ MT E M -i eg \\ ' e S E ,*A°) \\f E v e e \\ ' 1 J 



where the freedom of phase is fixed by Se,o(0) = 9 + k and S E ,ij>,o(0) = + k + if> as well as 

and z>f 



the continuity in v. Invariant measures v E , z>^ and v E for these actions and the probability 



measures p, p^ and p^ are defined by 

dv E (9)f(9) = [ du E (e)E v f(S E , v (e)mod7t) , /gC(R/ttZ), 



and similar formulas for z>^ and z>^. Due to a theorem of Furstenberg BLJ, the invariant 
measures exist and are unique whenever the associated Lyapunov exponent is positive. Let us 
note that one can easily verify that the invariant measure 0^ is simply given by the Lebesgue 
measure ^. Furthermore and v E do not depend on p (but v E does). 

Next let us write out a more explicit formula for the new Lyapunov exponent j^E). First 
of all, according to Furstenberg's formula |BL| IJSSj . 

loo(E) = S+ tV f di>g{e) log(||f|„ e ,||) , where e e = ( ^) ) ■ 

As is the Lebesgue measure, rotations are orthogonal and the integrand is 7r-periodic, one 

gets 

loo(E) = - E v J - \og[(ee\(MT E M^T(MT E M^)\ee)) . (4) 

Now (MT E M~ 1 )*(MT E M~ 1 ) = |1 + Pf| 2 is a positive matrix with eigenvalues A„ > 1 and 

1/ A„ where A^ = 1 + | + \J a> + with a = sin ^ fc ^ • As it can be diagonalized by an orthogonal 
transformation leaving the Lebesgue measure invariant, we deduce that 
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7oc(£) = \ ^J Q 2n ^ log(A,cos 2 W + i-sin 2 w) = IJ dp(v) log(^) ■ 
This formula shows immediately that ^^(E) > unless p = So (in which case A„ = Ao = 1). 

3 Result on the Lyapunov exponent 

Theorem Let E = —2 cos(fc) with either - rational or k satisfying the weak diophantine con- 
dition 

\l-e 2mk \ > ce~ eiml V mGZ, (5) 

for some c > and £' > 0. Then 

k satisfies (JSJ) , 

k = vr2 , 
depending only on p, 

\%(E)-^ foo (E)\ < ce-^\ . 

The result can be interpreted as follows: if the density of the impurities is small, then the 
incoming (Priifer) phase at the impurity is uniformly distrubuted for a sufficiently irrational 
rotation angle {i.e. quasi-momentum) because the sole invariant measure of an irrational rota- 
tion is the Lebesgue measure. For a rational rotation, the mixing is to lowest order in p perfect 
over the orbits of the rational rotation, which leads to the definition of the family (^^ J )(p (T ) e s 
and its distribution p q . As indicated above, the proof that this is the correct image is another 
simple application of modified Priifer phases and an oscillatory sum argument. 

Let us note that %(E) ^ ^(E); more detailed formulas for the difference are given below. 
As a result, one can expect a numerical curve of the energy dependence of the Lyapunov 
exponent at a given fixed low density to have spikes at energies corresponding to rational 
quasimomenta with small denominators. Moreover, the invariant measures v E and are 
not close to the Lebesgue measure, but has higher harmonics as is typical at Kappus-Wegner 
anomalies. Furthermore, let us add that at the band center E — the identity j(p, 0) = P7o(0) 
holds with no higher order correction terms and where % (0) is the center of band Lyapunov 
exponent of the usual Anderson model with distribution p. 

Finally, let us compare the above result with that obtained for a weak-coupling limit of 
the Anderson model [PF, JSS : the Lyapunov exponent grows quadratically in the coupling 
constant of the disordered potential, while it grows linearly in the density. The reason is 
easily understood if one thinks of the change of the coupling constant also rather as a change 
of the probability distribution on the space of matrices. At zero coupling, the measure is 
supported on one critical point (or more generally, on a commuting subset), and the weight in 
its neighborhood grows as the coupling constant grows. In ((TJ the weight may grow far from 
the critical point, and this leads to the faster increase of the Lyapunov exponent. 



P7 oo(£) + 0(p 2 



1(P,E) 



p%(E) + 0(?) 
where p and q are relatively prime. Furthermore, for £ 
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4 Proof 



For fixed energy E, configuration (v n ) n( z^ and (-0 n )neN, as well as an initial condition 9q, let us 
define iteratively the seqences 



n 



SE,v n {9n-l) , 9n = <SE,ip n ,v n (9n-l) = $E,v n (^n-l) + V'n ■ (6) 

When considered modulo tt, these are also called the modified Priifer phases. They can be 
efficiently used in order to calucate the Lyapunov exponent as well as the density of states. For 
the Lyapunov exponent, let us first note that one can make a basis change in (J2J) at the price 
of boundary terms vanishing at the limit, and furthermore, that according to |BL| A. III. 3. 4] it 
is possible to introduce an arbitrary initial vector, so that 



7 ( P ,£) = MiEbgNI^Ji^jr^lj 



(7) 



Now using the modified Priifer phases with initial condition 9 = 9, this can be developed into 
a telescopic sum: 



^ n=l 

1 N 

= p lnr^ - E E„ log ( MT^M^e^ 

n=l 

where in the second step we have evaluated the partial expectation over the last random variable 
v n by using the fact that for v n = the contribution vanishes. Next let us note that the function 
e 10 i — > E„ \og(\\MTy M _1 ee||) has an analytic extension to C\{0}, contains only even frequencies 
so that its Fourier series 



E„ log 



MT E M -l e& 



E 



dm, 6 
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has coefficients satisfying for any £ > a Cauchy estimate of the form 



Comparing with (jlj), we deduce 



Introducing now the oscillatory sums 



N 



a = ■y oc (E) . 



N 



UN) = E - , L(N) = E - J2 

n=l n=l 

the Lyapunov exponent now reads 

7(/),£) = /) V a m lim I m (N) . 



(8) 



(9) 
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Hence we need to calculate I m (N) perturbatively in p. Clearly Iq(N) = 1. Furthermore, 
integrating over the initial condition w.r.t. the invariant measure gives for all N 6 N 



dv E (9) I m (N) = / du E {9) e 



2im8 



Hence calculating I m (N) perturbatively also gives the harmonics of v E perturbatively (similar 
statements hold for I m (N), of course). Going back in history once, one gets 



N 



I m (N) = 4 E X) ((1 - P) e 2jmfc e 21 ™ 9 - 1 + p % i {^ mS ^ 9 ^fj 

n=l 

= (1 - p) e 2mk I m (N) + 0{p,N- x ). 
For k satisfying (J5J), one deduces 

\UN)\ < eM][ 0{ P ,N- x ) < ce^O(p,N-i). 

Replacing this and (JHJ) with £ > £' into (JHJ) concludes the proof in this case because only the 
term m = gives a contribution to order p. 

If now k = 7T-, the same argument implies 

Wr(JV) = O^N' 1 ) , VneZ, r = l,...,g-l. (10) 

Setting 

^^ e 2vmS E ,v(0)^ _ ^ fo( m ) g2t(m+Q6> E^e 2 " 7 ^*^) = g(m) e 2i(rn+l)9 

lez lei 
one deduces for the remaining cases 

I n? (iV) = (1 - p) I n? (iV) + ^(iV- 1 ) + p t q) (ln q+ i{N) + OiN- 1 )) . 

Due to (JTUj), this gives the following equations 

J n9 (iV) = I {n+r)q (N) + 0{{pN)-\p). 

They determine the invariant measure v E to lowest order in p. This shows, in particular, that 
v E is already to lowest order not given by the Lebesgue measure. We will not solve these 
equations, but rather show that the oscillatory sum I nq (N) satisfy the same equations, and 
hence, up to higher order corrections, the measure v E can be used instead of v E in order to 
calculate the Lyapunov exponent. Indeed, it follows directly from (0) and the definition of p q 
that 

tirn) r fS m ^ 

"i — c mmo dg=o ; 
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In particular, I m (N) = if mmodg 7^ 0. Thus we deduce 



L q (N) = + 0{N- 1 ) = i (n+r)q (N) + OiN- 1 ). 

Comparing the equations for I nq (N) and I nq {N) (which have a unique solution becaue the 
invariant measures are unique by Furstenberg's theorem), it follows that 

L q (N) = I nq (N) + 0(p, (pN)- 1 ) . 
Replacing this into (JDJ), one deduces 



j(p,E) = p ^ a m f dvf(6)e 2me + 0(p 2 ) 

= p J dv*(6)V v log (|| MT*M- l e e ||) + 0(p 2 ) . 

Now due to the orthogonality of rotations one may replace MT^M~ X by T£ v , and then the 
r.h.s. contains exactly the Furstenberg formula for 'jq(E) as claimed. The estimate comparing 
%(E) and 00(E) follows directly from the Cauchy estimate (JSJ). 

5 Result on the density of states 

Another ergodic quantity of interest is the integrated density of states, defined by 

N(p, E) = lim — E # {negative eigenvalues of the restriction of H^—E to £ 2 ({1, . . . , iV})}. 

Recall, in particular, that jV(0, E) = k if E = —2 cos(fc). Defining the mean phase shift of the 
impurities by 

(p{9) = E V (S E>V (6) - 6) , 
the low density expansion of the density of states reads as follows: 



Af(p,E) 



1 - p) k + pf g <p(9) + 0(p 2 ) k satisfies © 



{l-p)k + pfd0?{O)m + 0(p 2 ) k = n^ 

with p and q relatively prime. The proof of this is completely analogous to the above when the 
rotation number calucation (e.g. JSSJ for a proof) is applied: 



1 - 

N(p,E) = lim -E^fS^O-M 



n=l 



I N 

[l-p)k + p lim — E &(9n-i) 

iV— >oo 1\ ^— ' 



n=l 
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